Mapping disease regulatory circuits at cell-type resolution from single-cell multiomics data

Resolving chromatin remodeling-linked gene expression changes at cell type resolution is important for understanding disease states. We describe MAGICAL, a hierarchical Bayesian approach that leverages paired scRNA-seq and scATAC-seq data from different conditions to map disease-associated transcription factors, chromatin sites, and genes as regulatory circuits. By simultaneously modeling signal variation across cells and conditions in both omics data types, MAGICAL achieved high accuracy on circuit inference. We applied MAGICAL to study Staphylococcus aureus sepsis from peripheral blood mononuclear single-cell data that we generated from infected subjects with bloodstream infection and from uninfected controls. MAGICAL identified sepsis-associated regulatory circuits predominantly in CD14 monocytes, known to be activated by bacterial sepsis. We addressed the challenging problem of distinguishing host regulatory circuit responses to methicillin-resistant- (MRSA) and methicillin-susceptible Staphylococcus aureus (MSSA) infections. While differential expression analysis failed to show predictive value, MAGICAL identified epigenetic circuit biomarkers that distinguished MRSA from MSSA.

A robust, highly multiplexed mass spectrometry assay to identify SARS-CoV-2 variants

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants are characterized by differences in transmissibility and response to therapeutics. Therefore, discriminating among them is vital for surveillance, infection prevention, and patient care. While whole viral genome sequencing (WGS) is the "gold standard" for variant identification, molecular variant panels have become increasingly available. Most, however, are based on limited targets and have not undergone comprehensive evaluation. We assessed the diagnostic performance of the highly multiplexed Agena MassARRAY(R) SARS-CoV-2 Variant Panel v3 to identify variants in a diverse set of 391 SARS-CoV-2 clinical RNA specimens collected across our health systems in New York City, USA as well as in Bogota, Colombia (September 2, 2020 - March 2, 2022). We demonstrate almost perfect levels of interrater agreement between this assay and WGS for 9 of 11 variant calls ({kappa} [≥] 0.856) and 25 of 30 targets ({kappa} [≥] 0.820) tested on the panel. The assay had a high diagnostic sensitivity ([≥]93.67%) for contemporary variants (e.g., Iota, Alpha, Delta, Omicron [BA.1 sublineage]) and a high diagnostic specificity for all 11 variants ([≥]96.15%) and all 30 targets ([≥]94.34%) tested. Moreover, we highlight distinct target patterns that can be utilized to identify variants not yet defined on the panel including the Omicron BA.2 and other sublineages. These findings exemplify the power of highly multiplexed diagnostic panels to accurately call variants and the potential for target result signatures to elucidate new ones. ImportanceThe continued circulation of SARS-CoV-2 amidst limited surveillance efforts and inconsistent vaccination of populations has resulted in emergence of variants that uniquely impact public health systems. Thus, in conjunction with functional and clinical studies, continuous detection and identification are quintessential to inform diagnostic and public health measures. Furthermore, until WGS becomes more accessible in the clinical microbiology laboratory, the ideal assay for identifying variants must be robust, provide high resolution, and be adaptable to the evolving nature of viruses like SARS-CoV-2. Here, we highlight the diagnostic capabilities of a highly multiplexed commercial assay to identify diverse SARS-CoV-2 lineages that circulated at over September 2, 2020 - March 2, 2022 among patients seeking care at our health systems. This assay demonstrates variant-specific signatures of nucleotide/amino acid polymorphisms and underscores its utility for detection of contemporary and emerging SARS-CoV-2 variants of concern.

Intrahost evolution and forward transmission of a novel SARS-CoV-2 Omicron BA.1 subvariant

AO_SCPLOWBSTRACTC_SCPLOWPersistent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections have been reported in immune-compromised individuals and people undergoing immune-modulatory treatments. Although intrahost evolution has been documented, to our knowledge, no direct evidence of subsequent transmission and stepwise adaptation is available. Here we describe sequential persistent SARS-CoV-2 infections in three individuals that led to the emergence, forward transmission, and continued evolution of a new Omicron sublineage, BA.1.23, over an eight-month period. The initially transmitted BA.1.23 variant encoded seven additional amino acid substitutions within the spike protein (E96D, R346T, L455W, K458M, A484V, H681R, A688V), and displayed substantial resistance to neutralization by sera from boosted and/or Omicron BA.1-infected study participants. Subsequent continued BA.1.23 replication resulted in additional substitutions in the spike protein (S254F, N448S, F456L, M458K, F981L, S982L) as well as in five other virus proteins. Our findings demonstrate that the Omicron BA.1 lineage can diverge further from its already exceptionally mutated genome during persistent infection in more than one host, and also document ongoing transmission of these novel variants. There is an urgent need to implement strategies to prevent prolonged SARS-CoV-2 replication and to limit the spread of newly emerging, neutralization-resistant variants in vulnerable patients.

Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p=5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights. Author SummaryCOVID-19 clinical outcomes vary immensely, but a patients genetic make-up is an important determinant of how they will fare against the virus. While many genetic variants commonly found in the populations were previously found to be contributing to more severe disease by the COVID-19 Host Genetics Initiative, it isnt clear if more rare variants found in less individuals could also play a role. This is important because genetic variants with the largest impact on COVID-19 severity are expected to be rarely found in the population, and these rare variants require different technologies to be studies (usually whole-exome or whole-genome sequencing). Here, we combined sequencing results from 21 cohorts across 12 countries to perform a rare variant association study. In an analysis comprising 5,085 participants with severe COVID-19 and 571,737 controls, we found that the gene for toll-like receptor 7 (TLR7) on chromosome X was an important determinant of severe COVID-19. Importantly, despite being found on a sex chromosome, this observation was consistent across both sexes.

RT-PCR/MALDI-TOF diagnostic target performance reflects circulating SARS-CoV-2 variant diversity in New York City

As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to circulate, multiple variants of concern (VOC) have emerged. New variants pose challenges for diagnostic platforms since sequence diversity can alter primer/probe binding sites (PBS), causing false-negative results. The Agena MassARRAY(R) SARS-CoV-2 Panel utilizes reverse-transcription polymerase chain reaction and mass-spectrometry to detect five multiplex targets across N and ORF1ab genes. Herein, we utilize a dataset of 256 SARS-CoV-2-positive specimens collected between April 11, 2021-August 28, 2021 to evaluate target performance with paired sequencing data. During this timeframe, two targets in the N gene (N2, N3) were subject to the greatest sequence diversity. In specimens with N3 dropout, 69% harbored the Alpha-specific A28095U polymorphism that introduces a 3-mismatch to the N3 forward PBS and increases risk of target dropout relative to specimens with 28095A (relative risk (RR): 20.02; p<0.0001; 95% Confidence Interval (CI): 11.36-35.72). Furthermore, among specimens with N2 dropout, 90% harbored the Delta-specific G28916U polymorphism that creates a 3-mismatch to the N2 probe PBS and increases target dropout risk (RR: 11.92; p<0.0001; 95% CI: 8.17-14.06). These findings highlight the robust capability of Agena MassARRAY(R) SARS-CoV-2 Panel target results to reveal circulating virus diversity and underscore the power of multi-target design to capture VOC.

Acute COVID-19 gene-expression profiles show multiple etiologies of long-term sequelae

Summary ParagraphTwo years into the SARS-CoV-2 pandemic, the post-acute sequelae of infection are compounding the global health crisis. Often debilitating, these sequelae are clinically heterogeneous and of unknown molecular etiology. Here, a transcriptome-wide investigation of this new condition was performed in a large cohort of acutely infected patients followed clinically into the post-acute period. Gene expression signatures of post-acute sequelae were already present in whole blood during the acute phase of infection, with both innate and adaptive immune cells involved. Plasma cells stood out as driving at least two distinct clusters of sequelae, one largely dependent on circulating antibodies against the SARS-CoV-2 spike protein and the other antibody-independent. Altogether, multiple etiologies of post-acute sequelae were found concomitant with SARS-CoV-2 infection, directly linking the emergence of these sequelae with the host response to the virus.

Robust clinical detection of SARS-CoV-2 variants by RT-PCR/MALDI-TOF multi-target approach

The COVID-19 pandemic sparked rapid development of SARS-CoV-2 diagnostics. However, emerging variants pose the risk for target dropout and false-negative results secondary to primer/probe binding site (PBS) mismatches. The Agena MassARRAY(R) SARS-CoV-2 Panel combines RT-PCR and MALDI-TOF mass-spectrometry to probe for five targets across N and ORF1ab genes, which provides a robust platform to accommodate PBS mismatches in divergent viruses. Herein, we utilize a deidentified dataset of 1,262 SARS-CoV-2-positive specimens from Mount Sinai Health System (New York City) from December 2020 through April 2021 to evaluate target results and corresponding sequencing data. Overall, the level of PBS mismatches was greater in specimens with target dropout. Of specimens with N3 target dropout, 57% harbored an A28095T substitution that is highly-specific for the alpha (B.1.1.7) variant of concern. These data highlight the benefit of redundancy in target design and the potential for target performance to illuminate the dynamics of circulating SARS-CoV-2 variants.

Molecular Profiling of COVID-19 Autopsies Uncovers Novel Disease Mechanisms

BackgroundCurrent understanding of COVID-19 pathophysiology is limited by disease heterogeneity, complexity, and a paucity of studies evaluating patient tissues with advanced molecular tools. MethodsAutopsy tissues from two COVID-19 patients, one of whom died after a month-long hospitalization with multi-organ involvement while the other died after a few days of respiratory symptoms, were evaluated using multi-scale RNASeq methods (bulk, single-nuclei, and spatial RNASeq next-generation sequencing) to provide unprecedented molecular resolution of COVID-19 induced damage. FindingsComparison of infected/uninfected tissues revealed four major regulatory pathways. Effectors within these pathways could constitute novel therapeutic targets, including the complement receptor C3AR1, calcitonin-like receptor or decorin. Single-nuclei RNA sequencing of olfactory bulb and prefrontal cortex highlighted remarkable diversity of coronavirus receptors. Angiotensin I converting enzyme 2 was rarely expressed, while Basignin showed diffuse expression, and alanyl aminopeptidase was associated with vascular/mesenchymal cell types. Comparison of lung and lymph node tissues from patients with different symptomatology with Digital Spatial Profiling resulted in distinct molecular phenotypes. InterpretationCOVID-19 is a far more complex and heterogeneous disease than initially anticipated. Evaluation of COVID-19 rapid autopsy tissues with advanced molecular techniques can identify pathways and effectors at play in individual patients, measure the staggering diversity of receptors in specific brain areas and other well-defined tissue compartments at the single-cell level, and help dissect differences driving diverging clinical courses among patients. Extension of this approach to larger datasets will substantially advance the understanding of the mechanisms behind COVID-19 pathophysiology. FundingNo external funding was used in this study. Research in contextO_ST_ABSEvidence before this studyC_ST_ABSInformation regarding changes seen in COVID-19 has accumulated very rapidly over a short period of time. Studies often rely on examination of normal samples and model systems, or are limited to peripheral blood or small biopsies when dealing with tissues collected from patients infected with SARS-CoV-2. For that reason, autopsy studies have become an important source of insights into the pathophysiology of severe COVID-19 disease, highlighting the emerging role of hyperinflammatory and hypercoagulable syndromes. Studies of autopsy tissues, however, are usually limited to histopathologic and immunohistochemical evaluation. The next frontier in understanding COVID-19 mechanisms of disease will require generation of highly dimensional, patient-specific datasets that can help dissect this complex and heterogeneous disease. Added value of this studyOur work illustrates how high-resolution molecular and spatial profiling of COVID-19 patient tissues collected during rapid autopsies can serve as a hypothesis-generating tool to identify key mediators driving the pathophysiology of COVID-19 for diagnostic and therapeutic target testing. Here we employ bulk RNA sequencing to identify key regulators of COVID-19 and list specific mediators for further study as potential diagnostic and therapeutic targets. We use single-nuclei RNA sequencing to highlight the diversity and heterogeneity of coronavirus receptors within the brain, suggesting that it will be critical to expand the focus from ACE2 to include other receptors, such as BSG and ANPEP, and we perform digital spatial profiling of lung and lymph node tissue to compare two patients with different clinical courses and symptomatology. Implications of all the available evidenceCOVID-19 is a far more heterogeneous and complex disease than initially anticipated. Advanced molecular tools can help identify specific pathways and effectors driving the pathophysiology of COVID-19 and lead to novel biomarkers and therapeutic targets in a patient-specific manner. Larger studies representing the diversity of clinical presentations and pre-existing conditions will be needed to capture the full complexity of this disease.

Before the Surge: Molecular Evidence of SARS-CoV-2 in New York City Prior to the First Report

New York City (NYC) emerged as a coronavirus disease 2019 (COVID-19) epicenter in March 2020, but there is limited information regarding potentially unrecognized severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections before the first reported case. We utilized a sample pooling strategy to screen for SARS-CoV-2 RNA in de-identified, respiratory pathogen-negative nasopharyngeal specimens from 3,040 patients across our NYC health system who were evaluated for respiratory symptoms or influenza-like illness during the first 10 weeks of 2020. We obtained complete SARS-CoV-2 genome sequences from samples collected between late February and early March. Additionally, we detected SARS-CoV-2 RNA in pooled specimens collected in the week ending 25 January 2020, indicating that SARS-CoV-2 caused sporadic infections in NYC a full month before the first officially documented case.

Topoisomerase 1 inhibition therapy protects against SARS-CoV-2-induced inflammation and death in animal models.

The ongoing pandemic caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is currently affecting millions of lives worldwide. Large retrospective studies indicate that an elevated level of inflammatory cytokines and pro-inflammatory factors are associated with both increased disease severity and mortality. Here, using multidimensional epigenetic, transcriptional, in vitro and in vivo analyses, we report that Topoisomerase 1 (Top1) inhibition suppresses lethal inflammation induced by SARS-CoV-2. Therapeutic treatment with two doses of Topotecan (TPT), a FDA-approved Top1 inhibitor, suppresses infection-induced inflammation in hamsters. TPT treatment as late as four days post-infection reduces morbidity and rescues mortality in a transgenic mouse model. These results support the potential of Top1 inhibition as an effective host-directed therapy against severe SARS-CoV-2 infection. TPT and its derivatives are inexpensive clinical-grade inhibitors available in most countries. Clinical trials are needed to evaluate the efficacy of repurposing Top1 inhibitors for COVID-19 in humans.

Cytotoxic lymphocytes are dysregulated in multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and multiple organ involvement in individuals under 21 years following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. To identify genes, pathways and cell types driving MIS-C, we sequenced the blood transcriptomes of MIS-C cases, pediatric cases of coronavirus disease 2019, and healthy controls. We define a MIS-C transcriptional signature partially shared with the transcriptional response to SARS-CoV-2 infection and with the signature of Kawasaki disease, a clinically similar condition. By projecting the MIS-C signature onto a co-expression network, we identified disease gene modules and found genes downregulated in MIS-C clustered in a module enriched for the transcriptional signatures of exhausted CD8+ T-cells and CD56dimCD57+ NK cells. Bayesian network analyses revealed nine key regulators of this module, including TBX21, a central coordinator of exhausted CD8+ T-cell differentiation. Together, these findings suggest dysregulated cytotoxic lymphocyte response to SARS-Cov-2 infection in MIS-C.

Pathophysiology of SARS-CoV-2: targeting of endothelial cells renders a complex disease with thrombotic microangiopathy and aberrant immune response. The Mount Sinai COVID-19 autopsy experience

BACKGROUNDSevere Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) and its associated clinical syndrome COVID-19 are causing overwhelming morbidity and mortality around the globe, disproportionately affecting New York City. A comprehensive, integrative autopsy series that advances the mechanistic discussion surrounding this disease process is still lacking. METHODSAutopsies were performed at the Mount Sinai Hospital on 67 COVID-19 positive patients and data from the clinical records were obtained from the Mount Sinai Data Warehouse. The experimental design included a comprehensive microscopic examination carried out by a team of expert pathologists, along with transmission electron microscopy, immunohistochemistry, RNA in situ hybridization, as well as immunology and serology assays. RESULTSLaboratory results of our COVID-19 cohort show elevated inflammatory markers, abnormal coagulation values, and elevated cytokines IL-6, IL-8 and TNF. Autopsies revealed large pulmonary emboli in four cases. We report microthrombi in multiple organ systems including the brain, as well as conspicuous hemophagocytosis and a secondary hemophagocytic lymphohistiocytosis-like syndrome in many of our patients. We provide electron microscopic, immunofluorescent and immunohistochemical evidence of the presence of the virus and the ACE2 receptor in our samples. CONCLUSIONSWe report a comprehensive autopsy series of 67 COVID-19 positive patients revealing that this disease, so far conceptualized as a primarily respiratory viral illness, also causes endothelial dysfunction, a hypercoagulable state, and an imbalance of both the innate and adaptive immune responses. Novel findings reported here include an endothelial phenotype of ACE2 in selected organs, which correlates with clotting abnormalities and thrombotic microangiopathy, addressing the prominent coagulopathy and neuropsychiatric symptoms. Another original observation is that of macrophage activation syndrome, with hemophagocytosis and a hemophagocytic lymphohistiocytosis-like disorder, underlying the microangiopathy and excessive cytokine release. We discuss the involvement of critical regulatory pathways.

Introductions and early spread of SARS-CoV-2 in the New York City area

New York City (NYC) has emerged as one of the epicenters of the current SARS-CoV2 pandemic. To identify the early events underlying the rapid spread of the virus in the NYC metropolitan area, we sequenced the virus causing COVID19 in patients seeking care at the Mount Sinai Health System. Phylogenetic analysis of 84 distinct SARS-CoV2 genomes indicates multiple, independent but isolated introductions mainly from Europe and other parts of the United States. Moreover, we find evidence for community transmission of SARS-CoV2 as suggested by clusters of related viruses found in patients living in different neighborhoods of the city.